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Summary 

A quasistatic unilateral frictionless contact problem for a rigid axisymmetric indenter 
pressed into a homogeneous, linearly elastic and transversely isotropic elastic layer 
bonded to a homogeneous, linearly elastic and transversely isotropic half-space is 
considered. Using the general solution to the governing integral equation of the 
axisymmetric contact problem for an isotropic elastic half-space, we derive exact 
equations for the contact force and the contact radius, which are then approximated 
under the assumption that the contact radius is sufficiently small compared to the 
thickness of the elastic layer. An asymptotic analysis of the resulting non-linear 
algebraic problem corresponding to the fourth-order asymptotic model is performed. 
A special case of the indentation problem for a blunt punch of power-law profile is 
studied in detail. Approximate force-displacement relations are obtained in explicit 
form, which is most suited for development of indentation tests. 

1. Introduction 

In recent years, development of minimally invasive measurement techniques for determining 
the biomechanical properties of biological and artificial replacement tissues has become a 
problem of paramount importance in medicine. The mechanical properties of biomaterials 
can be evaluated by means of simple indentation techniques ([T]), which have been proved 
useful for both identification of mechanical properties of a soft biological tissue, like articular 
cartilage ((2]), and assessing its viability during arthroscopy ([3]). Indentation tests are 
performed with various indenter geometries, and among the most frequently used types 
of indenters are axisymmetric (in particular, cylindrical, spherical, and conical) indenters. 

Development of an indentation test for the articular cartilage layer requires solving 
the corresponding contact problem for an elastic layer (see Fig. [iji). The widely used 
mathematical model for indentation tests of articular cartilage developed by Hayes et al. 
([2]) is based on the analytical solution (|4| of the frictionless contact problem for an elastic 
layer bonded to a rigid base. The Hayes et al. model allows to take into account the 
thickness effect, but, at the same time, it neglects the substrate effect. In order to take 
into account the influence of the substrate deformation, we consider the frictionless contact 
problem for an elastic layer bonded to an elastic half-space (see Fig.!!]}). 
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Fig. 1 a) Schematics of indentation test for the articular cartilage layer attached to subchondral 
bone; b) Schematics of the indentation problem for an elastic layer bonded to an elastic half-space. 



In the present paper, we consider the quasistatic unilateral contact problem for a 
rigid axisymmetric indenter pressed into a homogeneous, linearly elastic and transversely 
isotropic elastic layer bonded to a homogeneous, linearly elastic and transversely isotropic 
half-space. We focus on deriving asymptotic approximations for the force-displacement 
relation (that is the relation between the contact force P and the indenter displacement w), 
which constitutes the basis for developing depth-sensing indentation tests. In the isotropic 
case, the substrate effect in the flat-ended and spherical indentation was studied recently in 
([5]) based on the asymptotic solutions obtained in J6l [7]) . The main novelty of the present 
paper consists in considering the case of an arbitrary axisymmetric indenter. As a special 
case, we study the indentation problem for a blunt punch of power-law profile. 

The main difficulty for an analytical analysis of the axisymmetric unilateral contact 
problem is associated with the fact that the radius a of the circular contact area u> is 
not known in advance and must be determined as a part of the solution. Apparently the 
first asymptotic analysis of the axisymmetric contact problem for an arbitrary axisymmetric 
indenter dates back to England, whose work ([8]) was published in 1962. As a special case, 
England (|8]) considered the indentation problem for a hemispherically ended punch and 
obtained asymptotic expansions in terms of powers of a small parameter /i = R/h, where 
R is the curvature radius of the hemispherical end, h is the thickness of the elastic layer. 

Our asymptotic analysis is based on a small parameter e = a/h defined as the ratio of 
the contact radius to the layer thickness. It should be noted that in contrast to the method 
originally developed by England (8), we construct asymptotic expansions in terms of an 
unknown quantity in the same way as in the asymptotic method of large A developed by 
Vorovich et al. ([6]). However, our approach differs from previous studies ((6l[9| in two ways. 
First, we do not construct an asymptotic expansion for the contact pressure, but for the aim 
of deriving the force-displacement relation, following (|10l Hip , we derive exact equations 
for the contact force and the contact radius and operate with integral characteristics of 
the contact pressure. Second, we perform asymptotic analysis of the resulting non-linear 
algebraic problem obtaining the approximate force-displacement relation in explicit form, 
which is most suited for development of indentation tests. 
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2. Indentation problem formulation 

We consider the unilateral frictionless contact problem for an elastic layer bonded to an 
elastic half-space, which is made of a different material. Let the composite elastic medium 
occupying the half-space x 3 > be indented by an arbitrary frictionless axisymmetric 
indenter. Assuming that the a^-axis coincides with the axis of symmetry of the indenter, 
we will describe its surface by the equation 

x 3 = $>(xi,x 2 ), (2.1) 

where <&(xi,X2) is the shape function whose dependence on the Cartesian coordinates X\ 
and xi comes through a polar radius r = \j x\ + x\. Without loss of generality, we may 
assume that <3>(0, 0) = 0. 

Under the assumption that the indenter shape function §{x\,xi) is convex and the 
indenter is loaded with a normal force, P, which is directed along the positive a^-axis, 
contact between the indenter and the surface of the elastic medium will be established over 
a circular region, w, say, of radius a. For a blunt indenter, the contact area uj is variable 
and depends on the contact force P. 

Further, let Gs{x\, x 2 , 0) be the surface influence function for the elastic medium that 
gives the vertical surface displacement of the elastic layer under a unit point force applied at 
the origin of coordinates and directed along a^-axis. Then, the unilateral contact problem 
under consideration, which involves an a priori unknown contact area w, reduces to that of 
solving the following integral equation: 



G 3 (xx - ?/i,a:2 - 2/2,0)p(y) dy = w - $(xi,x 2 ). (2.2) 
Here, p(xi,X2) is the contact pressure, w is the indenter displacement. 



It should be emphasized that Eq. (2.2 1 implicitly contains an a priori unknown radius a 
of the contact area w, which is to be determined as a part of solution from the following 
conditions: 

p(y) > for < |y| < a, p(y) = for |y| = a. (2.3) 

Here, |y| = \Jy\ + y| is a polar radius. 

By the equilibrium equation, the external contact force P is related to the contact pressure 
density p{x\, x-i) as follows: 

' r p(y)dy = P (2.4) 



Due to the axisymmetry of the elastic medium, the surface influence function can be 
represented as 

G3(yi,y 2 ,0) = -L( * (2 _ 5) 



2tt6» V |y| h \ h 

where h is the layer thickness, 9 is an elastic constant related to the first layer, F(t) is a 
dimensionless function. 

It can be shown (j6]) that the following integral representation holds: 

oo 

F{t) =/"[!- L(u)] J (ut) du. (2.6) 
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Here, Jo (a;) is the Bessel function of the first kind, the kernel function L(u) depends both 
on the material properties of the elastic layer and the elastic half-space. 

According to (jfij, there is a neighborhood of zero on which the function (2.6) is expanded 
into an absolutely convergent power series 

oo 

F{t) = a m t 2m (2.7) 

m=0 

with the coefficients 

oo 



[(2m) 



1112 



J [1 - L{u)]u 2m du. (2.8) 



Thus, the indentation problem is to find the force-displacement relationship, P(w). In 
the present paper, we employ an asymptotic approach to derive approximate analytical 
formulas for P(w) in the small-contact limit, that is under the assumption that the contact 
radius a is small compared with the layer thickness h. 

3. Equation for the radius of the contact area 

In the axisymmetric case, Eq. (2.2 1 can be reduced to a simpler form (see, e.g., (|6l l9f) 



where K(s) is the complete elliptic integral of the first kind, J 7 (a, r) is a dimensionless 
function given by 

/•oo 

T{<t,t)= I [1- L(u)]Ja(au)J (Tu)du. (3.2) 
Jo 

In the neighborhood of the origin, the function J- (a, r) has the following absolutely 
convergent expansion (|9|): 

oo oo 

J( ff ,T) = ^^^, (3.3) 

i=0 3=0 

(to') 2 



where a m are given by (2.8 1. 

Now, introducing an auxiliary function v(r) by the formula 



v(r) = 9(w - $(r)) + - / T (|, £j p{ P ) P dp (3.4) 



and making use of the general solution of the axisymmetric contact problem for an elastic 
isotropic half-space (|12U13l [14|) . we represent the solution to Eq. (3.1) in the form 



ASYMPTOTIC ANALYSIS OF THE SUBSTRATE EFFECT 



■5 



7rV(r) = v(0) + r 



v (s) ds 
\Jr 2 — s 2 



(3.6) 



From the boundary condition (2.3 1 of vanishing the contact pressure at the edge contour 

(3.7) 



of the contact area, it is necessarily follows that 

V(o) = 



Hence, in view of (3.7), formula (3.5) reduces to 

, , i r v>( s )<is 

P{r) 



7T J \/s 2 — r 2 



(3.8) 



Now, according to Eq. (3.4) and the assumption $(0) = 0, we will have 



2?r 



(3.9) 



(3.10) 



The substitution of the expression (3.2) into Eqs. (3.9), (3.10) and then into Eq. (3.6) 
yields the following formulas (|10|) : 



V(r) = 2dw - 20r 



&(s)ds 2 



V V 2 — s 2 h 



So (~, ^ P{p)pdp, 



So{cr, t) = / [l — L(u)\ Jo(au) cos rudu. 



(3.H) 



(3.12) 



<$>>{p)dp 
y/a 2 - p 2 Oh 



Thus, Eq. (3.7) for determining the contact radius a takes the form 

So(Z~)p(p)pdp. (3.13) 

It is clear that if an approximation for the contact pressure p(r) is known, the second 
term of the right-hand side of Eq. (3.13 1 will give a correction to the corresponding equation 
of the Galin — Sneddon theory of axisymmetric elastic contact (|15l 1161) . 

4. Equation for the contact force 



(2.4), we will have 



Substituting the expression (3.8) for the contact pressure into the equilibrium equation 

a 

P = 2 J V(s) ds. 



(4.1) 
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Now, the substitution of (3.11) for V(s) into Eq. (4.1) yields 



a a 

P = A6wa-A6 J &{p)^a 2 - p 2 dp + A J S 2 



p a 
h' h 



p(p)pdp, 



where 



S2(cr,a) = j [l — L(it)] Jq((tu) du 



(4.2) 



(4.3) 



Finally, taking into account Eq. (3.13), we eliminate the variable w in Eq. (4.2) as follows: 

(4.4) 



P = 40 



&{p)p 2 dp 



4 / S- 



(^l)p(p)pdp, 



Here we introduced the notation 

oo 

S 2 {(J,a)= / [l - L(u)] Jq(ctu) 



sin au 



a cos au du. 



(4.5) 



Again, if we know have an approximation for the contact pressure p(r), then the second 
term of the right-hand side of Eq. (4.4) will produce a correction to the corresponding 
equation of the Galin — Sneddon theory of axisymmetric elastic contact (|15l 1 1 6 [) . 

5. The fourth-order asymptotic model 

According to Eqs. (3.12) and (4.5), the following series expansions hold: 



oo oo 



i=0 j=0 



(5.1) 



(0) 2 2 ^(m!) 2 

« = (z!) 2 (2j)! ftm ( m = z + ^)' 



OO oo 



»=o j=i 



(5.2) 



, (2) 2 2 ^j(m\) 2 
« " (i!) 2 (2j + l)! 



a m (m = i+j). 



Now, keeping terms of the infinite series (5.1 1 and (5.2) which contain only the coefficients 
ao and oi, we get 

^(cr, a) ~ a + ai(er 2 + 2a 2 ), (5.3) 



S 2 (<7, a) ~ --aia 3 . 



(5.4) 
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The substitution of the asymptotic expressions (5.3) and (5.4) into Eqs. (3.13) and (4.4) 
results in the following approximate relations: 



1 f P ( 

v/a 2 - p 2 tor9h\ V 



a Q + 2ai 



/i 2 



-Pi 



«i 
/i 2 



P~ 40 



V a2 ~ P 



1^- VF* 



(5.5) 



(5.6) 



Here, P and P 2 are, respectively, the contact force and the polar moment of inertia of the 
contact pressure, i.e., 



P = 2nj P(p)pdp, 


a 

p 2 =2^y p 2 p{p)pd P . 



(5.7) 



Substituting the expression (3.8) for the contact pressure into Eq. (5.7), we obtain 

a 

p 2=^J V(s)s 2 ds. 



(5- 



Now, the substitution of (3.11) for V(s) into Eq. (5.8) yields 



w - / 



\/a 2 - p 2 



\ 3 



Here we introduced the notation 

oo 

5 4 (<t, a) = / [l - L(u)l Jo(ctu) 



sin aw 2 / sin aw 



aw (ait) 2 V ait 



/smait \ . 

cos au } du 

V ait / 



(5.10) 



In view of the fact that the quantity P 2 enters Eq. (5.5) with the coefficient a\, Eq. (5.9) 
may be reduced to the following one being in the same range of accuracy with Eqs. (5.5) 

(2p 2 - a 2 ) , \ , 4a 3 a 



and (5.6): 



P-2 



w- / $(p) 



\/a 2 ~ p- 



rpdp 



P. 



Note that in writing down Eq. (5.11 1, it is assumed that 

S 4 (cr,a) ~ ^a . 



(5.11) 



(5.12) 
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Thus, the three relations (5.5), (5.6), and (5.11) form the resulting algebraic problem 
with respect to the four variables w, P, P2, and a. 



Eliminating the variable P2 from Eq. (5.5) by the substitution (5.8), we obtain 

a a 

4ai a 3 ' 



w 1 + 



3tt h 3 



<$>'{p)dp Aai_ 
» irh 3 



P 
2n8h 



a Q + 2ai 



,, Mp 2 -a 2 ) 

Hp) - r-n — k p d P 



\faP r - p 2 



h 2 



4aoai a 
3tt P 



On the other hand, Eq. (5.6) yields 



(5.13) 



Pi 1 



8ai a 3 
3tT/i3 



40 



&(p)p 2 dp 
\/a 2 - P 2 



(5.14) 



With the same asymptotic accuracy as that of Eqs. (5.5), (5.6), and (5.11), we obtain 



from Eq. (5.14) the following approximate relation: 

8fl i fl3 \ „„ 1 &{p)P 2d P 



P 



1 - 



3tt h 3 



y/a 2 - p 2 



(5.15) 



Now, substituting the expression (5.15) into Eq. (5.13), we arrive after some asymptotic 



simplifications at the following expansion: 

<&'(p)dp 4ai 



w 



1 



4ai a 3 
37 h 3 



\Ja 2 - p 2 



ixh 3 



Mp 2 ~a 2 ) 

Hp) r-n — k P d P 



2 

nil 



a Q + 2a 1 



a 2 8a a 1 a 3 \ f <&'{p)p 2 dp 



h 2 



3tt h 3 



(5.16) 



Thus, Eqs. (5.15) and (5.16) define the force displacement relation in a parametric way. 

We note that the fourth-order asymptotics of the axisymmetric unilateral contact problem 
for an elastic layer was previously derived by Alexandrov and Pozharskii ([9| in a somewhat 
different way (see, ([9]), Ch. 1.3, Eqs. (25) and (26)). It can be shown that the two solutions 
are asymptotically equivalent up to terms of order 0(e 5 ) as e — > 0, where e = a/h. 



6. The case of a blunt indenter 

Let us consider indentation of the elastic medium by a blunt indenter with the shape function 

$(r) = Ar\ (6.1) 

where 1 < A is an arbitrary real number, A is a constant having the dimension [L 1 ~ A ] with 
L being the dimension of length. 
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Taking into account the integral 

1 



/A+1 '" 2A ~ 2A rf A ) 2 r ( A)- 



vT^F (A + l) V2 



we reduce Eqs. (5.15) and (5.16) to the following ones, respectively: 



P ~ 6>AFi(A)a A+1 1-e 



,8ai\ 
3W : 



w~ AF 2 (X)a x < l-e 



2Xa 3 8A(2A + 5)ai 



16Aaoai 

'tt(A + 1) " 3tt(A + 3)(A + 1) ' " 3vr 2 (A + 1) 
Here we introduced the notation 

-Ax 2 



\2oA 

f 1 (x) = A^r 



(A + l) V2 
A\2 



r(A)- 1 , 



F 2 (A)=A2 A - 2 r(-) r(A)" 1 , 



(6.2) 
(6.3) 

(6.4) 

(6.5) 
(6.6) 



Following (17) and using the asymptotic method based on Lagrange's formula for solving 



algebraic equations (18J, we invert Eq. (6.3) to find the contact radius as a function of the 
indcnter displacement in following form 



- ~ mil + Byuo + Bnw 2, + B^w 6 + B 4 zu 4 
h v 



(6.7) 



with the coefficients given by 



Si 



2a 



tt(A + 1)' 



Bo 



(A + 3) /2a y 

2(A+l) 2 l 7T J 



Bo 



(A + 4)(A + 2) /2a \3 (2A + 5) 8oi 
3(A + 1) 3 V 7T ) + (A + l)(A + 3) 3tt 



5 4 



(2A + 5)(A + 5)(3A + 5)(A + 3) /2a \ 4 (A 2 + 11A + 22) 16a a 



24(A + l) 4 (A + 3) 
Here we introduced the notation 



\ tt J (A + l) 2 (A + 3) 3tt 2 



1 / w \i 
h\AF 2 (X)) 



Now, substituting the expansion (6.7) into Eq. J6.2L we obtain 



where 



P ~ 6A-^F 3 (X)w^(l + C x w + C 2 vj 2 + C 3 m 3 + C 4 m 4 

'X- 



A-l A + 2 



F 3 (A) = X — 2—(X + l)- L T 



r(A): 



(6.8) 

(6.9) 
(6.10) 
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Ci 



C, = 



2oo „ _ (2A + 3) (2a y 

TT ' 2 2(A + 1) V 7T J ' 

(A + 2)(3A + 4) /2a \3 (A + 2) 801 
3(A + l) 2 V TT / + (A + 3) 3tt ' 



(2A + 5)(3A + 5)(4A + 5) /2a \ 4 (2A + 5)(A + 2) 16a a 
w — . TVo + 7^ . -i \ / a . „n ^r^- 



24(A + 1) 3 



V TT / (A + l)(A + 3) 3tt 2 



(6.11) 



In view of (6.8), Eq. (6.9) represents the load-displacement relationship for a blunt 
indenter with the shape function (6.1 1. Note that Eqs. (6.7) and (6.9) generalize the second- 



order asymptotic model developed in (|19(1 . 
6.1 Example. Paraboloid indenter 



In the case A = 2 and A — 1/(2R), Eq. (6.1 ) takes the form 

,.2 



$(r) 



r 
2R' 



(6.12) 



where R is the radius of curvature of the indenters surface at its vertex. 
According to Eqs. (6.4), (6.5), (6.8), and (6.10), we will have 

4^2 



Fi(2) 



16 



*s(2) 



Bi = 



c x = 



2a 
3tt : 



2a 



Bo = 



10a 2 

97T2 



B, = 



_ 320a 3 , 
63 -^7^ 



64ag 8ai 
27tt 3 + W' 

32ai 



B 4 = 



154ag 
27tt 4 



256aoai 
45tt 2 : 



15tt 



C 4 = 



286ao 64a a! 



9tt 4 



+ 



5tt 2 



The obtained results are in agreement with those in (|19p . 
6.2 Example. Cone indenter 

For a cone indenter, we have A — 1 and A = tan 7, whereas Eq. (6.1 ) takes the form 

$(?') = r tan 7, 

where 7 is the angle between the contact surface and the side surface of the cone. 



According to Eqs. (6.4), (6.5), (6.8), and (6.10), now we have 

ir 4 
, ^3(1) = - 



*!(!)= TT, F 2 (l) 



B x = 



r; 



p _ 2«g 



B* = 



7 a 1 
3P 



5 4 = ^ 



34aoai 



2w cot 7 
nh ' 
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C X = 



2an 



r _ 5ao 

<^2 — 5": 



c 3 = 



14a 3 , 2ai 



C 4 = 



42ap 14a ai 



Substituting the above coefficients into Eqs. (6.7) and (6.9), we arrive at the following 
approximate relations: 



2 cot 7 f 2an w 8a 2 . , w 2 
M 1 + ^ cot 7 - + -f cot 2 7T ^ 



7T 7T* 

/40a 3 , ^ 56ai 
1, ! 



3tt 4 



COt d 7TT + 



/224a^ 544a ai\ 



h 3 \ 7T 6 



3tt 6 J 



cot 7 



/* 4 ' 



(6.13) 



P 



40 cot 7 9 f „ 4a w 20a 2 , , 

'ii; 2 <^ 1 + -^cot 7 - + — ^cot 2 7 

7r y tt z h 7r 4 



UT 



112a 3 , 



16aA 



C0t + 



/672a^ 224a ai 



/l 3 \ 7T 8 



cot 7 



h* 



(6.14) 



Equations (6.13) and ( 6.14[ ) generalize the second-order asymptotic model developed in 
()19|) . We note the following misprints in (|19]k missing factors cot 7 and cot 2 7 in Eq. (31) 
and an extra factor § in Eq. (32). 



7. The case of a hemispherically-ended indenter 

For a hemispherical indenter, we will have 



(7.1) 



where i? is the curvature radius of the hemispherical end. 
In what follows, we make use of the identities 



&(p)dp 1, fl + a 
= ^ In 



2 _ p 2 2 



1 — a 



p 



...» I ol 1 o. /I + a 
2 =i? i-2 + 4 (1 + a)ln 



(2p 2 -a 2 ) 

i/a 2 - p 2 



pdp = R 4 



12 



1~ 4,, /1 + a 



1 - a 



; (l-a 4 )ln 



Here we introduced the notation 



a 

5' 



Thus, Eqs. (|5.15| and (|5.16|) can be represented as follows: 

'1 + a 



P ~ 6>i? 2 (1 + a 2 ) In 



1 - a 



2a 



o o8ai 



(7.2) 



(7.3) 
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w 
R 



In 



1 + a 
I- a 



-^((l + a 2 )ln 



1 + a 
I- a 



3°i 



Gtt 

4 4a ai 3 



2a(7a 2 + 3) - (7a 4 + 6a 2 + 3) In 
1 + a 



- 2a 

1 + a 
1 - a 



3tt 2 



a d [ (1 + a 2 )ln 



1 - a 



2a 



Here we introduced the notation 



R 
h' 



(7.4) 



(7.5) 



Now, following England ([8]), we construct the approximate solution for the contact radius 
in the form of a power expansion with respect to the small parameter [i as 



R 



a + fioei + /-t 2 «2 + M 3q: 3 + A* «4- 



(7.6) 



The coefficients of the expansion (7.6) are determined according to Eq. (7.3). By a 
perturbation method, we get 



a 3 



a% =0, a 2 — 0, 04 = 0, 
4fl! « 2 (1 - a 2 )[(l + a 2 ) ln(±±f*) - 2a ] 



3tt (i_ a g )ln (i±m) +2ao 
while ao is a unique root of the equation 



P 
OR* 



(l + «8)la(i±^) 



2a . 



In view of (7.7), it is readily seen that the expansion (7.6) simplifies to 

a 



R 



a + fi 6 a 3 + 0(jj, b ). 



(7.7) 
(7.8) 

(7.9) 

(7.10) 



The substitution of (7.10) into Eq. (7.4) yields the following result after expanding up to 
the fourth order in /1: 



w 
R 



In 



1 + a Q 
1 - a a 



~ M g((l + a 2 )lu 



1 + a 



2a, 



V?M (3a 4 , + 2a 2 + 3) In 

D7T 



1 + ap 
l-a 



1 - ao, 

-6a (l + «o) ) TOW'). (7.11.) 



Formulas (7.6) -(7.11) are in complete agreement with the corresponding formulas from 

2a 



{[8]) , provided the following relations hold: 



fci 



1 + k\ 1 + ^2/71-72 



1 



1 



K a = 



K, = 



4ai 



Here, k\ 1 k%, 71, 72 are elastic constants used in (8), Kq and K\ are the corresponding 
asymptotic constants. 
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8. Indentation scaling factor for a blunt indenter 

Let us rewrite Eq. (6.9) in the form 

P = 6A-iF 3 (\)w^ 1 f x (w), 



introducing the indentation scaling factor f\(w). According to (6.9), we have 
f\(zu) — 1 + C\w + C^nj 2 + Ojnj 3 + C4CC7 4 , 



where the coefficients C n are given by (6.11). 

Observe that the factor fxfa) depends on the dimensionless indentation parameter w 
defined by Eq. (6.8). In view of Eqs. (6.3) and (6.8), we have 

2a _ 2 ( x ~ l ) ( 2a o \ 2 

7T(A + 1) £ 2(A+1) 2 V 7T J 

_ 3 A A - 1 )( 2A - 1 )/2ao\ 3 



6(A + 1) 3 



8ai(2A + 5) 
3tt(A + 3)(A + 1) 



(A-l)(2A-l)(3A-l) /2a \ 4 16a ai (A 2 - A - 8) 



24(A + 1) 



(v) 



3^ 2 (A + l) 2 (A + 3) 



(8.3) 



Now, substituting the expansion (8.3) into (8.2), we express the indentation factor in 

(8.4) 



terms of the relative contact radius by introducing a new notation 

«a(e) = /a Me)). 



Thus, from (|8.2|) and (|8.3|), it follows that 

«a(e) 



1 + £ 2oo +£2 (2^1)^200^2 



2(A + 1) V 7T J 
3 /(2A+l)(3A + l) /2a \ 3 8oi(A- 2) 



6(A + 1) 2 V 7T J ' 3tt(A + 3) 
(2A+l)(3A + l)(4A + l)^2a ^ 4 , 16a ai (2A 2 + 4A - 1) 

7T 



24(A+ l) 3 

Substituting the value A = 2 into Eq. (|8.5|), we readily get 



3tt 2 (A + l)(A + 3) 



(8.5) 



«2(e) 



2a 2 10 a o 

7T 37T 2 

4 /70aQ 16aoai\ 



, / 140a 3 , 32ai 
1 27tt 3 + 157 



V 9tt 4 



3^ r 



In order to check the asymptotic formula ( 8.5 ) , we pass to the limit as A 
the following result: 



(8.6) 

-00, obtaining 



Kqo (e) 



1 + e 



2a 



2ao 

7T 



4 32aoai\ 



8ai 
3^ 



3^ )■ 



(8.7) 



The case A = +00 corresponds to a cylindrical indenter. The asymptotic expansion ( |8.7[ ) 
is in complete agreement with the fourth-order asymptotics derived previously in (j6l I20[) . 
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9. Displacement-force relationship for a blunt indenter 



According to the notation ( |6.8[ ), the force-displacement relation (6.9 1 can be rewritten as 
P ~ 6AF 1 (\)h x+1 w x+1 {\ + dm + C 2 vo 2 + C 3 m 3 + C 4 zu 4 ), (9.1) 



(9.2) 
(9.3) 

(9.4) 



where the factor Fi(X) is given by ( |6.4[ ). 
Inverting the relation (9.1), we obtain 

w ~ P(l + DiP + D 2 P 2 + D 3 P 3 + DiP 4 ) , 

where we introduced the notation 

P= (9AF 1 (\))~^h~ 1 P^ r[ . 



Now, in view of (6.8), from Eq. (|9.2|) it follows that 
P 



w ~ A*+ 



h( ^ 2Aa ° 
\9F 3 (X)J \ ir(\ + l 



r(A 



-P 



8A(A + 2)ai 
3tt(A + 1)(A + 3) j 



Formula (9.4) represents the sought-for relation between the indenter displacement w and 
the contact force P. 



between the contact radius a and the contact force P in the simple form 



Substituting now the expansion (9.2) into Eq. (6.7), we derive the corresponding relation 

(9.5) 



a 
h 



P 1 



3tt(A + 1) 



P 



Observe that Eq. (9.5) can be directly obtained from Eq. (6.2). 
10. Incremental indentation stiffness for a blunt indenter 

During the depth-sensing indentation, the incremental indentation stiffness can be evaluated 
according to the relation 

dP (10.D 

da 



dP 
da 

dw ~ dw ' 



Substituting the expressions (6.2) and (6.3) into Eq. (10.1), we arrive at the approximate 
relation 



dP 

dw 



49a- 



1 



8ai (A+4) 3 
3tt (A+l) fc 



(10.2) 



' 1 _ 2oo _ _ 8ai (2A+5) 3 , 16a ai (A+4) ^4 ' 
1 7T b 3tt (A+l) fc I" 3tt 2 (A+l) fc 

Note that in deriving formula ( 10. 2[ ), the following relation was taken into account (sec 

(A + l) Fx (A) 



Eqs. (RU} and (fgTEt): 



= 4. 



A F 2 (A) 

Now, expanding the right-hand side of Eq. ( |10.2[ ) into a power series in e, we obtain 



dP „ { 2a Q „ / 2a 

— ~ AOall + e— +e 2 ( — 

dw 7T V 7T 



3 f {^o\ 3 , 801 



V 7T / 



3tt 



32aoai 
3tt 2 



(10.3) 
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It is interesting to observe that the coefficients of the asymptotic expansion (10.3) do not 
depend on the parameter A, which describes the indenter shape. 

On the other hand, the incremental indentation stiffness for a blunt indenter can be 



evaluated by differentiating the force-displacement relationship (9.1| with respect to the 
indenter displacement. Thus, by taking into account the formula 



we get the relation 

dP 
dw 



dm zu 1 x 
dw~ = XF^XyAh 1 



) (A + ^ F 3 {X)F 2 (X)^hzu\\ + wd (A 



A + l 



2 r (A + 3) 

ZU C-2" 



A + l 



rSr (A + 4) 

ZU G3- 



A + l 



ir (A + 5) 

ZU O4 



A+l 



(10.4) 



where the coefficients C„ are given by (|6.11 1. 
Now, substituting the expression 



terms of powers of the parameter e, we again arrive at Eq. (10.3) by noting that 



3 ) for zu into Eq. ( 10.4 ) and expanding the result in 



(A + l) 
A 



F 3 (X)F 2 (Xp =4. 



Finally, by comparing Eqs. (10.3) and 
rewritten in the following form: 

dP 

dw 



7), it is readily seen that Eq. (10.3) can be 



46»a/i oc (£). 



(10.5) 



Here, Kqo(£) is the indentation scaling factor for a cylindrical indenter. 



Observe that Eq. ( 10.5 1 can be regarded as a generalization of Barber's theorem ()21l) for 



the incremental indentation stiffness established in the case of an isotropic elastic half-space. 



Note also that in context of indentation testing, Eq. (10.5) should be rewritten as 




^oo(£): 



(10.6) 



where A is the current area of contact. 

In contrast to the Bulychev — Alckhin — Shorshorov (BASh) relation established in (|22l 
1231 1241) in the case of an isotropic elastic half-space, Eq. (10.6) contains the extra factor 
«oo(e), which accounts for both the thickness and substrate effects. Note that in the case 
of isotropy, we have 6 = E^/[2(l - (i/ 1 )) 
and Poisson's ratio of the layer material. 



')], where E^ and are Young's modulus 



11. Conclusion 

An asymptotic analysis of the substrate effect has been performed for deriving approximate 
solutions of the quasistatic indentation problem for an arbitrary frictionless indenter. 
Explicit asymptotic formulas are derived in the framework of the fourth-order asymptotic 
model, which has been previously shown to be effective for describing the thickness effect for 
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a spherical indenter on an elastic layer. The results obtained can be applied for development 
of indentation tests. 
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Appendix: Kernel function for a transversely isotropic layer bonded to a 
transversely isotropic half-space 

The constitutive relationship for a transversely isotropic material referred to the Cartesian 
coordinates (xi, £2, £3) with the Ox\Xi plane coinciding with the plane of elastic symmetry 
can be written in the matrix form as follows (25): 



( au ) 




" A n 


A12 


A13 













{ 


£11 


\ 


C22 




Ayi 


An 

















£22 




CT33 






A 13 


^33 















£33 




013 













2A44 












£13 




C23 
















2A44 









£23 




V ^12 / 

























£12 


/ 



For a transversely isotropic material, only five independent elastic constants are needed to 
describe its deformational behavior. The elastic moduli An, A12, A13, A33, and A44 can 
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be expressed in terms of the engineering elastic constants as follows (26): 



El 



An = 



e_ 12 

E' u 



(1 + ^)1 



2E v n 

E< " 



E v 



A 12 = 



E' u 



(1 + ^)1 



Ay-i = 



Ev' 



A iZ = 



E'{1 - u) 



Aaa = G'. 



Here, E and E' are Young's moduli in the plane of transverse isotropy and in the direction 
normal to it, respectively, v and v' are Poisson's ratios characterizing the lateral strain 
response in the plane of transverse isotropy to a stress acting parallel or normal to it, 
respectively, G' is the shear modulus in planes normal to the plane of transverse isotropy. 
Note also that A n - A 12 = E/(l + v). 

In the case of a transversely isotropic elastic layer bonded to a transversely isotropic 
clastic half-space, in accordance with the known solution (]27|) (see Eqs. (102) and (103)), 
we will have 

1 



L(u) = l + g(exp(- U / 7 f ) ),exp(- U /7^ ) )), 

where 

2M{x uyi ) 
Q{xi,yx) = — r . 

According to Eq. (101) (27), the functions M{x\ : yi) and N{x±,yi) are given by 

mi \ ( (i) (i) ^ (l) 2 , (i) 2 

M(xi,yi) = -(g\' 012 - g 2 a 2 i)Xiy 1 - g\ 'anx 1 + g\ ' a 22 y x 

- (012021 - aiia 2 2)x\yl, 

N{xi,yi) = 1 + 2(g[ 1) a 12 - g^a^xxyx + (g[ X) + g^^ianxl - a 22 yf] 

+ (012021 - ana 22 )xlyf. 
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According to Eqs. (51) -(54) (|27l) . the coefficients an, a\2, 021, and a 2 2 are defined as 



an = 1 + 



27i 



(i) 



Z 



-H 2 



(2) / (: 
7i - ("M j 7 2 



( 2 )\ 2 „, ( 2 ) 



gH 



2( 7 f 3 -m^ 7 ^) 



(2)„,(2)^ 



+ (m^-l)(m( 2 )-l) 7 ^ 



„(!) 



(2) (2) 

7i - 72 



012 



27 



(i) 



«21 



27i 



(i) 



^2^(!) nTTlnW (2)W (1) -.N 



,(2) 



(2) (2) 



022 = 1 



27: 



(l) 



Z 



(2) / (2)\2 (2) 

7i ^ (™i ) 72 



(mi i) -l)(mf ) -l) 7 l 



.,(!) 



1 



(2) (2) 

7i - 72 



2mP( 7i °-m< 2) 7 f) 

+ g 2 (™[ 1} f 



Here the following notation is used (see Eqs. (48), (59), and (60) (|27|) ): 
Z = ^{( 7 f)- 7 «)( ff f)-(mf)) 2 5 f)} 



(2) (2) 

7i - 72 
2f 7 f ) - m 



(1) (2) i\/ (1) (1) , (2) (2)\ 

i\ - 1) (m\ - 1) (7} + 7} ^ ; J 



f) 7 f)( 



(i) (i) (1) 
7i - m\ ; 7 2 



(1) / (i)\ 2 (l) 
7i ~ \jn\ ') 7 2 



Si 



(i) 



7i 



(i) 



(1) (i)' 
7i - 72 



(1) 

^2 (a) 

(i) (1)' »i 
7i - 72 



(2) 

Tj (2) 
(2) (2) ' ^2 

7i - 72 



72 



(2) 



(2) (2) ' 

7i - 72 



-H = : — —s -, g ~ 



ffWfm^-l) 



(2) (2) ' 

7i - 72 



We note two misprints in formula (59) (|27p for Z: one closing bracket is missing and the 
product gH is misprinted as H g . 

The layer material is characterized by the elastic constants H^\ mj[ , 7} , 72 , an d 
so on, while the material of the elastic half-space is characterized by the elastic constants 
H^ 2 \ m[ , 7 j 2 ' ) , 7 2 2 \ and so on. So, we have 



H (n) = 



)= ( 7 ^+ 72 -y 
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and 



(n) _ A ll \ll ) 



A 



(n) 
13 



A 



A 

W 

44 



(n) 



(n) _ A ll yTg ) 



A 



(n) 
13 



A 



A 

44 



00 



Finally, the dimensionless parameters 7} , 73 are denned as the roots of the equation 



4 jWiW _ ^2 

7 ^11 ^44 7 



A u ^33 



.4 



3 ^13 



2A 



^33 ^44 



(11.1) 



In the case of isotropy, we have 

L(u) = l + Q(u) ) Q(u) = -2e~ 2u 



di + d 2 e~ 2u 
1 + d 3 e- 2u + d 2 e,- 4u '' 



2(1 -(vny) 

where 2?W and i/ 1 ) are Young's modulus and Poisson's ratio of the layer material, while 
the coefficients di, d 2 , and d 3 are given by formulas (105) -(108) (27). We note that in 
this particular case, Fabrikant's solution obtained in (|27p is in agreement with Burmister's 
solution (1281). 



